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Abstract. We report the measurement of the leptonic (e“''e“) decay channel of the 
0(1020) meson in Au-|-Au collisions at = 200 GeV by the STAR experiment. 

The transverse momentum (pt) spectrum is measured for 0.1 < px Y 2.5 GeV/c 
at mid-rapidity (|j/| < 1). We obtain the px-integrated 0(1020) mass = 1017.7 ± 

0.8(stat.)±0.9(sys.) MeV/c^ and width P^ = 8.0±2.5(stat.)±2.3(sys.) MeV/c^, which 
are within 1.5 cr and 1.1 u of the vacuum values, respectively. No significant difference 
is observed in the measured px spectrum, dN/dy, or (px) of the 0(1020) meson 
between the e+e“ and hadronic {K~^K~) decay channels as measured by the same 
experiment. The experimental results are compared to a theoretical model including 
medium-modified 0(1020). 


PACS numbers: 25.75.-q, 25.75.Dw, ll.30.Rd, 13.20.-v 
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1. Introduction 

Experiments at Relativistic Heavy Ion Collider (RHIC) have found a variety of 
indications for the deconhnement phase transition and formation of Quark-Gluon 
Plasma (QGP) pQ. Chirality, a fundamental symmetry of QCD which is spontaneously 
broken in vacuum, is predicted to be restored within an energy density range similar 
to the deconhnement phase transition [2]. The high energy density and temperature 
achieved in relativistic heavy-ion collisions provide a unique environment to study the 
chiral phase transition, which is of crucial importance in the Standard Model. The 
properties of hadronic resonances, which are excitations of a (partially) chirally restored 
medium, are measured to study chiral symmetry restoration in hot and dense matter 
because short lifetimes allow a signihcant fraction of those produced to decay in-medium. 

Among resonances, the 0(1020) meson has several attractive features which can be 
used to study the two phase transitions. The 0(1020) meson, which is an ss bound 
state, is an ideal probe for the strangeness production in the medium [3]. As a vector 
meson, it can decay not only into kaons, but also into dileptons via a virtual photon 
state. Absence of strong interactions with hnal state hadrons in the dilepton channel 
permits reconstruction of 0(1020) decays from throughout the collision evolution. The 
mass, width, and branching ratios of 0(1020) are predicted to be modihed in a nuclear 
medium HEIEIE]. Since the mass of 0(1020) is just above the mass of two kaons, a 
decrease in the 0(1020) mass may lead to a change in the relative production rates of 
0(1020) —)■ and 0(1020) — K^K~ due to phase space limitations in the hadronic 

channel. Even though the lifetime of 0(1020) in vacuum (44 fm/c |8]) is longer than the 
estimated lifetime of the QGP (4-10 fm/c [H]), interactions with the medium may lead 
to an increase of the width, effectively shortening the lifetime. 

2. Analysis 

In this paper, we report the measurement of 0(1020) production via the e’''e“ decay 
channel in Au-|-Au collisions at = 200 GeV with the STAR experiment [TU] 

at RHIG. The results presented in this paper are consistent with a recent STAR 
publication within systematical and statistical uncertainties, although used analysis 
methods differ m- The detector subsystems used in this analysis are the Time 
Projection Ghamber (TPG) [12] and the Time-of-Flight (TOE) [T3] . 

2.1. Event Selection 

Events are selected with collision vertices within ± 30 cm of the center of the TPG along 
the beam axis ( 14 ), and within 2 cm transverse of the beam axis to exclude interactions 
with the beam pipe. To suppress background from overlapping collisions in the TPG 
and ensure that the TOE start time is calculated for the triggered collision, the Vertex 
Position Detector (VPD) [IT], which provides the TOE start time, is used to select 
events for which \VzypD — 14,tpc| < 3 cm. After an additional selection of the 0-80% 
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most central events (minimum bias for this analysis) based on the multiplicity of charged 
particles at mid-rapidity [T5|, about 250 million events are analyzed in this study. 


2.2. Track Selection 

In order to reduce the contamination by daughter particles from weak decays (which 
typically have a decay length of several cm) and the tracks from 7 conversions in the 
material of the beam pipe, we require the distance of closest approach to the primary 
collision vertex of tracks to be less than 1.1 cm (where the resolution is typically 0 . 2 - 0 .3 
cm). Tracks are required to have at least 22 points in the TPC. The ratio of number of 
measured to possible TPC points in the TPC must be greater than 0.52 to avoid split 
tracks. The transverse momenta (pt) of the tracks have to be larger than 0.18 GeV/c 
to reach the TOP detector, which is mounted outside of the TPC. The momenta of the 
electrons (including positrons if not specihed) are required to be less than 2 GeV/c due 
to the large identihcation contamination by charged hadrons above 2 GeV/c. 



A(3-V(3-^ (TOF) 


Figure 1. (Color online) The distribution of ncr(dE/dx) vs. for electrons, 

pions, and muons within a low momentum range. The red ellipse shows the cut to 
select electron candidates. 


2.3. Electron Identification 

To achieve best separation of electrons from other particle species, the velocity 
measurement fi from the TOF and the ionization energy loss (dE/dx) measurement 
from the TPG are used simultaneously. We apply particle identihcation (PID) cuts 
in a two-dimensional (2D) space: n(T(dE/dx) = {X — with X = In (dE/dx), 

which is a Gaussian distribution with mean px and sigma cxx, and j = 
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{l^TOF — /3g /fdrpQp, which is the deviation from electron expectation. First, we select 

electrons reqniring the distance from the mean of the electron distribntion in the 2D 
space within 3 ct (the red ellipse in Fig. [^, where a for both dE/dx and TOF depend on 
the momentnm and psendo-rapidity {rj) of the measured particles. In order to further 
suppress pion contamination, we apply an additional particle selection criteria based on 
the 2D probability [16] that a particle is an electron, which is defined as 


_ X Pdf(X; 0^) X Pdf(F; 0^) 

^ PDF(X;0)f) X Pdf(F;00)’ 


where Pdf is a probability density function and O^xy parameter sets for particle i. 

is the fraction of yield of particle i (^- iV* = 1). All the parameters and the fraction 
of yields are obtained by fitting particle distributions in the 2D space. The summation 
in the denominator includes only electrons and pions for simplification in this analysis. 
All Pdfs are Gaussian functions except the A/3~^//3~^ distribution for pions, for which 
the Student’s t function im is used to better describe the tails of the distribution. In 
this analysis, a probability larger than 60% is required for a particle to be identified 
as an electron. The purity of the selected electrons is estimated to be about 95% in 
0.18 GeV/c < Pt < 2.0 GeV/c by studying distributions of all relevant particle species 
such as pions, kaons and protons. 


2.4- (t)(1020) Signal Extraction 

The invariant mass of 0(1020) is reconstructed from all combinations of identified 
electron-positron pairs which have opening angles wider than 30 degrees in the same 
event (unlike-sign signal), shown as the black curve in the inset of Fig. The mixed- 
event technique [IHl HH] , which pairs electrons with positrons from 20 different events, 
is adopted to estimate the combinatorial background. In order to preserve the event 
characteristics in the mixed-event technique, e"*" and e~ pairs are selected only from the 
same event class, as defined by 14 position and event plane angle [20] with ten bins each, 
and eight centrality bins. After normalizing the mixed-event background to the unlike- 
sign signal distribution in the invariant mass regions of 0.91-0.95 GeV/c^ and 1.11-1.15 
GeV/c^ (orange areas in the inset of Fig. |^, it is subtracted bin by bin from the unlike- 
sign signal distribution to extract the 0(1020) signals. The resulting raw distribution 
for pairs with \y\ < 1 and px = 0.1 — 2.5 GeV/c is shown in Fig. [^as black points. A 
combined fit to the raw distribution, shown as the black curve in Fig. [^ is performed 
using a Voigt function to describe the signal and a quadratic polynomial to represent 
the residual background. The raw 0(1020) yield is the sum of the counts in each bin 
above the residual background fit for the range — 30 < Mee < -|- 30 MeV/c^ 

(which typically accounts for 85% of the total yield) and the integral of the peak fit 
function outside that signal range. We obtain a signal significance, defined as the raw 
signal yield over statistical error, of 15.9. To measure the px spectrum we divide the 
integrated signal into seven px bins with a significance of about 6-7 in each px bin. 
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2.5. Efficiency Correction for the px Spectrum 

The pt spectrum is corrected for the TPC track reconstruction efficiency, the matching 
efficiency between TPC tracks and TOP hits, the detector acceptance, the track quality 
cut efficiency, and the PID cut efficiency [16]. The track reconstruction efficiency and 
the TPC acceptance are obtained by embedding 0(1020) —>- 6 ^ + 6 “ tracks from Monte 
Carlo simulation into real Au+Au collision data at the detector response level and 
reconstructing them along with real tracks, such that intrinsic detector resolutions 
and inefficiencies are taken into account. Along with the intrinsic resolutions from the 
simulation, we further smear the simulated electrons’ momenta by cXp.^, = 0.51% px based 
on studies of the invariant mass line shape of J/ 0 , where we do not expect any medium 
modihcation to the width. The ratio of number of 0(1020) mesons after quality cuts to 
number of simulated 0 ( 1020 ) mesons is taken as the efficiency of the reconstruction with 
quality cuts. A sample of tracks from 7 conversions in the beam pipe (identified with 
89% purity) is used to estimate the TOP matching efficiency times acceptance. The PID 
cut efficiency is calculated directly from the two-dimensional probability distributions 
mentioned above [16] . The overall efficiency times acceptance is about 9% for the lowest 
px bin and about 4% for the highest px bin. 

2.6. Mass and Width Extraction 

To extract mass and width from invariant mass signals we use a Voigt function IZII. 
which is a Breit-Wigner function convolved with a Gaussian function accounting for the 
detector resolution, defined as 

/ OO 

dM'G{M;M',a)BW{M';Mo,T), (2) 

■OO 

where G{M-, M', a) is a Gaussian function with mean M' and sigma a, and 
BW{M'] Mq,T) is a non-relativistic Breit-Wigner function with mass Mq and width 
P. To resolve degeneracy between the resolution parameter a and the 0(1020) width, 
a simulation with no intrinsic P broadening is used to determine the value of a (i.e. P 
can be fixed to the vacuum value [22]). It is found that 1.16 MeV/c^ of additional px 
smearing beyond the intrinsic simulated resolution is necessary for a good £t, resulting in 
a value of 8.1 MeV/c^ for a, and further adding 0.04 MeV/c^ to the extracted mass Mq. 

This width extraction method is different from the recent STAR publication where 
the reported width includes the detector resolution. Taking the difference into account, 
the two results on the 0 ( 1020 ) width are consistent within statistical and systematical 
errors. 

2. 7. Systematical Errors Estimation 

The systematical errors are estimated by considering uncertainties due to: (1) 
mismatches between the real and simulated track distributions which lead to 
uncertainties in the track reconstruction efficiency, (2) the PID selection criteria, mainly 
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due to uncertainty in the 2D probability function, (3) the TOF matching efficiency, 
and (4) variations due to the choice of the signal extraction method, such as the 
fit functions for the signal (non-relativistic & relativistic Breit-Wigner function), the 
residual background function (first & second order polynomials), and the normalization 
ranges for the mixed-event background. These systematical uncertainties on yield, mass, 
and width are averaged over all pt bins and summarized in Tab. [Tj The systematical 
uncertainty due to the additional px smearing is estimated as 0.2 MeV/c (2.5%) in 
width, which is insignificant within the statistical and other systematical errors of our 
measurement and not added in the results. 


Category 

Yield 

Mass 

Width 

TPG track quality cut 

9.3% 

0.083% 

10 .0% 

Particle identification (PID) 

4.6% 

0.043% 

5.9% 

TOF matching efficiency 

8.3% 

0.0 

0.0 

Inv. mass signal extraction 

3.8% 

0.007% 

3.4% 

Additional px smearing 

0.0 

0.0 

2.5% 


Table 1. Systematical uncertainties averaged over px bins, shown as percentages of 
the measured values. 


3. Results 

3.1. Mass and Decay Width 

The px-integrated 0(1020) invariant mass distribution and fit result are shown in 
Fig. 1^ The extracted fit parameters are = 1017.7 ± 0.8(stat.) ± 0.9(sys.) MeV/c^ 
( 1 . 5(7 away from the Particle Data Group (PDG) value of 1019.5 ± 0.02 MeV/c^ [H]) 
and P</, = 8.0 ± 2.5(stat.) ± 2.3(sys.) MeV/c^ ( 1 . 1(7 away from the PDG value of 
4.27 ± 0.03 MeV/c^). We can extract an upper limit of 13.6 MeV/c^ on the width 
and a lower limit of 1015.7 MeV/c^ on the mass with a 95% confidence level. 

Figure [^shows the extracted mass and width of 0(1020) in each px bin from fitting 
the data. The mass deviates by 1.4 a and 1.6 a for px = 1.6 — 2 GeV/c and 2 — 2.5 
GeV/c respectively from the simulation, which exhibits the expected vacuum values in 
the detector. The deviations in width are 1.6 a and 2.1a for px = 0.4 — 0.7 GeV/c and 
1.6 — 2 GeV/c respectively. Most deviations are less than 2 a, but the measured masses 
are generally lower and the widths are systematically higher than vacuum values. We 
do not observe any px dependence for mass and width within their errors. 

3.2. Comparison to Hadronic Decay Channel 

The corrected 0(1020) —)■ invariant yields at rapidity \y\ < 1 for minimum 

bias Au-|-Au collisions at = 200 GeV are presented in Fig. along with 
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xIO^ 



Figure 2. (Color online) Black points are the px-integrated ((i(1020) —>■ e+e” signal 
from data after subtracting mixed-event background. The black curve is a fit to the 
signal with the Voigt function plus a quadratic polynomial, which represents residual 
background (orange long-dashed curve), and the gray band is a 95% confidence band 
of the fit. Red points and the short-dashed curve are the simulated ((>(1020) signal with 
the additional momentum smearing and the fit function to the simulation, respectively, 
with the residual background from fitting the data added. (Inset) Invariant mass 
distribution of the e“'"e“ pairs in the same event shown as a black curve, and 
normalized mixed-event background shown as a red curve. The orange areas show 
the normalization regions for the mixed-event background. Fit parameter units are 
GeV/c2. 


0(1020) — K^K~ results measured by STAR |23] as well as mx-exponential fits [23] . 
Both spectra are divided by the corresponding branching ratios to facilitate comparison. 
The inclusive 0(1020) yield per unit rapidity dN/dy is calculated by counting signals 
within the measured transverse momentum range and integrating the mx-exponential £t 
function for the rest. The mean transverse momentum (px) is obtained in the same way. 
The contributions from the unmeasured px regions, estimated via the extrapolation of 
the fitted function, are 3.7% and 8.3% for dN/dy and (px), respectively. We obtain 
dN/dy = 2.91 ± 0.16(stat.) ± 0.17(sys.) and (px) = 1-03 ± 0.06(stat.) ± 0.06(sys.) 
GeV/c. The values of dN/dy and (px) are in agreement with the values for measurements 
from the hadronic decay channel {dN/dy = 2.68 ± 0.15, (px) = 0.962 ± 0.015 GeV/c) 
within the errors (0.83cr and 0.79cr respectively). In the bottom panel of Fig. the 
relative differences between the two spectra, (We — Nkk)/Nkk, are shown. The Nkk 
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Figure 3. (Color online) Mass and width extracted from fitting as a function of p^. 
Dashed lines show the PDG ^(1020) mass and width values. Fit results to the data are 
shown as blue stars, whose statistical and systematical errors are represented by bars 
and filled boxes respectively. Red points are results of fits to the simulation with the 
additional momentum smearing (not shown for widths, as the fits use the fixed PDG 
value). The systematical uncertainties on the width from the simulation are included 
in those of the real data. 


are calculated by integrating the fit function of the hadronic decay channel result (dashed 
curve) in each px bin. The px spectrum of 0(1020) measured in the decay channel 
is in agreement with that from the hadronic decay within the statistical and systematical 
uncertainties. Similar results have been observed in the dimuon decay channel at SPS 
energies 

3.3. Comparison to a Theoretical Model 

We compare the measurements with a model calculation [25], which is based on 
an effective chiral hadronic Lagrangian with coupling constants determined from the 
hadronic vacuum properties. This model includes an in-medium 0 spectral function 
folded over the phase-space hreball evolution up to freeze-out [251 EZ]- One of two 
components in the model is a “cocktail” contribution, which comes from the 0(1020) 
decays after elastic hadronic interactions stop (kinetic freeze-out). The other component 
is from 0(1020) decays in a thermal medium (thermal radiation) with a modihed 0(1020) 
spectral function. Its broadening implies back reactions regenerating 0(1020) in the heat 
bath due to detailed balance. In this model, only 0(1020) decay after hadronization is 
included. Theoretically there can be an additional contribution, which is missing in this 
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Figure 4. (Color online) The corrected ^(1020) —>■ e+e" (blue star) and 0(1020) —?► 
K~ (red circle) [23] yields for the minimum bias events in Au + An collisions at 
= 200 GeV. The vertical bars are statistical errors and the boxes are systematical 
uncertainties. Systematical errors for the K~^K~ channel are fixed to be 10% for all 
Pt- The dashed and solid curves are mx-exponential fit functions to the hadronic and 
leptonic decay channel results, respectively. The gray band represents errors of the fit 
function to the leptonic result. Data points are placed at the mean px in each px bin 
estimated from the fit functions. 

model, from medium-modified 0(1020) decays in the QGP phase, if the 0(1020) bound 
state can exist at temperatures above the critical temperature of the QGP. For this 
comparison, the STAR detector conditions (resolution, acceptance and efficiency) are 
applied to the model predictions. 

The resulting invariant mass shape from the model is shown in Fig. In the model 
calculation, the in-medium radiation component shows a mass shift of 7 MeV/c^ and 
a width broadening of ~55 MeV/c^. Fitting the measured invariant mass distribution 
using the theoretical line shape, where the fraction of the in-medium contribution is hxed 
to be 10%, gives a of 1.40. However, by letting the fraction as a free parameter, 

the £t to data suggests a 29 ± 15% medium contribution with a y^/nd/ = 1.34. 


4. Summary 

In summary, we have presented the measurements of 0(1020) yield, mass, and width 
in the decay channel at mid-rapidity {\y\ < 1) in 0-80% minimum bias Au-|-Au 
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Figure 5. (Color online) Model fit to the (^(1020) —>■ e+e” px-integrated (over 
Pt = 0.1 — 2.5 GeV/c) invariant mass data, showing individual components of the 
model. The black curve is the Voigt fit function shown in Fig. 


collisions at = 200 GeV. The measurements of 0(1020) are consistent with its 

vacuum mass within 1.5 a and width within 1.1 a with no signihcant px dependence. 
No signihcant difference is observed in the measured 0(1020) px spectrum, dN/dy, 
or (px) between the leptonic (e’''e“) and hadronic {K~^K~) decay channels. The 
measured invariant mass distribution is compared to a model calculation including a 10% 
contribution from medium-modihed 0(1020) and a 90% contribution from late decays in 
the vacuum. Our data demonstrate modest sensitivity to the medium-modihed 0(1020) 
component, but insufficient to observe dehnitively the small integrated ehects predicted 
by this model. 
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